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Abstract. We prove a partial result concerning the long-standing problem 
on limit periodicity of the Jacobi matrix associated with the balanced measure 
on the Julia set of an expending polynomial. Besides this, connections of the 
problem with the Faybusovich-Gekhtman flow and many other objects (the 
Hilbert transform, the Schwarz derivative, the Ruelle and Laplace operators) 
that, we sure, are of independent interest, are discussed. 

1. Introduction 

In 80's the following interesting phenomena was discovered: the spectral measure 
of an almost periodic Jacobi matrix can be singular continuous (supported on a 
Cantor type set of the zero Lebesgue measure). The effect was studied from both 
sides — from coefficient sequences to spectral data ^ and from spectral data 
to Jacobi matrices. 

The second, usually more elegant, approach produced the following example ^j, 
0. Let T{z) = z'^ - C. For C > 2 the Jufia set of T is a real Cantor type 
set, \E\ = 0. Denote by fi the balanced measure on E, ^{T~^{F)) = ^{F) for all 
F CE. Let 

(1) J 



qo pi 
pi qi P2 



be the Jacobi matrix associated to the given measure. Note that to construct 
J : P{'L+) one uses the three term recurrent relation for polynomials 

orthonormal in L^^^ 

(2) \Pk{\) ^ PkPk-i{\) + qkPk{X) + Pk+iPk+i{X). 

of course Pk ^ where {|fc)} is the standard basis in /^(Z+). 
Then the given matrix satisfies the renormahzation equation: 

V*T{J)V = J, 

where V\k) = \2k). In fact, this is a system of nonlinear equations for p„'s (q„ — 
in this case), due to which at least for C > 3 one gets inductively that 

b2"i+m -Pm| < e„, for aU l,m\ e„ -> (n -> cx)). 

That is the sequence {pn\ and, by definition the matrix itself, is limit periodic. It 
looks very natural to conjecture that if only T is an arbitrary expanding polynomial 
in the sense of Complex Dynamics [7] then its balanced measure produces a limit 
periodic Jacobi matrix. Several research groups attacked this problem (in full 
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generality) but failed. Even the case of the quadratic polynomial with C > 2 
is still open. 

Recall some properties of Jacobi matrices. Let J be a Jacobi matrix, J* = J, 
acting in or P(Z+). Under the assumption pk 7^ the vector |0) of the standard 
basis is cyclic for J. The resolvent function is a function of the form 

(3) r{z) = {0\{J-z)-'\0). 

It has positive imaginary part in the upper half plane and hence possesses the 
representation 

(4) r{z)= /-^ = (l|(A-z)-i|l)^. . 



A 

where A- is the operator multiplication by the independent variable in and 1 
is the function that equals one identically. Formulas lj2Jl and Q give one to one 
correspondence between triples {L'^^n ^'^ 1} ^^^^ |0)} or {C'^, J, |0)}, re- 

spectively, in the finite dimensional case. To restore J starting from the nonnegative 
measure a one uses 

Our first object is the following 

Conjecture 1.1. Let T{z) be an expanding polynomial of degree d with a real Julia 
se<£;, T-i : ^ Define J = J{x) by 

(5) (o|(z- j(x))-i|o) = |^|^, xeK,e]. 

Respectively Jn{x) is associated with an iteration Tn — T°", degr„ = d„. Then for 
every e there exists n such that 

(6) J„(0)|| <e. 

Note that eigenvalues of J„(a;) and J„(0) are close, so the non trivial part deals 
with eigenvectors. 

Let us explain how this conjecture is related to the general one. If /i is the 
balanced measure on E, then the resolvent of J = J{lJ-) satisfies to the following 
Renormalization Equation 

(7) V*{z ~ J)-W = {T{z) ~ J)-'T'{z)/d, 

where V\k) — \kd). Let us include J into a chain {Jn{t)}te[o,i] defined by 

v:{z - j„(i))~ v„ = (r„(z) - tj)'^T:^{z)/dn 

(compare the last equation with ISJ). Then the main goal is to show that 

J„(0)|| <e forn>no, 

since it would imply immediately that J(/i) is limit periodic. Thus to prove Con- 
jecture ^31 is a good model problem on the way to prove limit periodicity of J(/i). 

The following approach looks very natural: to get ((HJ we have to estimate J'(x). 
The given derivative has a special representation 

(8) ^^FiJ) + [G,J] 

dx 

with F{J) = {T'{J)}~^ . It is a certain flow on Jacobi matrices that in a sense is 
dual to the well-known Toda flow. We call it FG flow (see Sect. 2). The first 
term at the right hand side in is small due to the characteristic property of 
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expanding polynomials: |T,'(a;)| > Ac''\ x E E, with A > 0, c > 1. It appears that 
the estimation we get for G is not enough to state that the commutator [G, J] is 
sufficiently small (Proposition [^21 • However on this way we found quite designing 
formulas and connections with so many objects (the Hilbert transform, the Schwarz 
derivative, the Ruelle and Laplace operators) that, we sure, they are of independent 
interest. 

In the framework of this approach, initiated in fP, we managed to prove the 
following theorem that partially confirms the main hypothesis. 

Theorem 1.2. Let J be the Jacobi matrix associated with iterations of an expanding 
polynomial T . Then for every e there exists n such that 

(9) \Pk+sdl -Pk\ < e, \qk+sdl - Qkl < e, 

for all s > and k = 1,2, d„. 

Note that actually our goal is to prove © when fc = l,2,...,d^. A proof of the 
theorem is given in Sect. 3. 

Acknowledgment. We wish to thank Misha Shapiro who called our attention to 
the results of 



2. FG FLOW 

2.1. Definition. Let J : C'' ^ C^. Consider the resolvent function 

(10) (o|(--^r^|o) = E7^- 

fc=i 

Under the Toda flow the spectrum is stable Afc — Const but masses vary with time 
CTfc = crfc(i). In FG flow case Afc = Afc(t) but CTfe — Const. Moreover, in our case (O 
time is x, ~ l/d and T(Afc(a;)) = x. Recall T{z) is an expanding polynomial of 
degree d with a real Julia set E, E C [-^,^1, : -> [~^,^]- 

We want to get a differential equation on J. Let 05 be a unitary matrix such 
that 

m = Q5A, 
where A = diag{Afc}. Since we can choose 

Ai(x) < A2(a;) < ... < Xd{x) 
that holds for all x, 23 essentially is well defined. We put 

■ Po(Ai) ... Po(Ad) ' 

Pd-i{X,) ... Pd-i{Xd)_ 

where Pk{z) is the orthonormal polynomial. 
We differentiate J with respect to x 

j = <BA<B"i + SBAS"^ - SAS^^SBS"^ ^ F + GJ - JG, 

where F := SA<B-\ G := ^^-^ . By the definition F = /(J) with /(Afc) = Afc. 
Thus F = T'(J)"^ The next step is to determine G. 
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Note some evident facts. G is skew-symmetric and (0|G — 0, so G|0) = 0. Also 
it is easy to show, say by induction, that 



r ^nr-^F + Gr 

ax 



Finally, since Pk{J)\Q) = \k) and 



-^Pfe(J) - |-Pfc(J) - FPliJ) + GPk{J) - Pk{J)G 
ax ox 



we get 
(11) 



d_ 

dx 



Pk{jm^FPi,{jm+GPk{jm. 



Let G_|_ be a lower triangle matrix with zeros on the main diagonal such that 
G = G+ - G;. Then ^ implies 

(12) G+|fc) = G+Pfe(J)|0) ^ -(FF^(J)|0))f . 

Here /i^-* means that in a vector ft, = {^j}j=o we have to replace all coordinates 
hj, < j < k, by zeros. 

Let us rewrite (I12f) in other words. Define an operator D by 

D\k)=DPk{J)\0) :=P{.(J)|0). 

Then 

G+ - -(PI?)+. 

It is easy to check using the functional representation in L^^ that 

(13) DJ-JD = I-\{pdPdy){Pd-i\, 

where PdPd{^) is defined by jSjl. Note that PdPd{z) — in L^^, that is it has the 
same roots {Afe(x)} as T{z) — x. Thus PdPd{z) ~ C{T{z) — x) and {pdPd)'{z) = 
GT'iz). 

Definition 2.1. FG flow is given by a differential equation of the form 

(14) j ^ F + GJ - JG 

with F — f{J) and G = G+ — G^, where G+ = —{FD)^ and D is an (upper 
triangle) matrix such that commutant \D^J] equals the unity matrix up to a one 
dimensional perturbation "H". 

2.2. {FD) as a Hilbert transform. 

Lemma 2.2. The matrix of the operator {FD) with respect to the basis of eigen- 
vectors of J has the form 



(15) 



1 T"(Ai) 

2 T'(Ai) 



Ai —Ad 



1 T"(Arf) 

2 T'(Ad) . 



T'(Ai) 



T'(Ad). 



Proof. Let us evaluate _D in the basis of eigenvectors of J. In this basis 
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As we know {pdPd)' {^k) — CT'{Xk). Taking into account that now J is diagonal 
we conclude that the diagonal entries of DJ — JD are zeros. Therefore Pd~i{Xj) — 
d/{CT'{Xj)}. Thus the right hand side of lO is of the form 

[T'iXi] 



I - 



1 1 



[T'{\d)\ 

Referring again to a diagonal form of J we solve and get 

1 r'(AO 



(16) 



D,, 



'■^ Xi~ XjT'{Xj) 
To find diagonal entries Da we have to use D\Q) = 0. Since 



1 



we get 



Note that 



Therefore 



-Di, = r'(A,)5;— i 



T'(Afe) Afe - A, 



1 

2-^ Ti(\, \ \, - 



k=l 



T'{Xk)Xk-z T{z)~x' 



(17) 



— — lim { — 
T'(AO .-A 



1 1 



1 



lim 



r'(A.) A, - z T{z) - X 

T(z)-x 1 -, A 



T{z)-x {T'{X,)f 
Thus (^Hll and finish the proof. 
2.3. Trace of {FD)*{FD). 

Lemma 2.3. Let L2 be a Ruelle operator of the form 
(18) L,g{x) ^\Y.{j^) 



Tv=x 



id let S{T) be the Schwarz derivative of T , S(T) = 



rp'n 

— 



itr{(FD)*(FD)}^-iL2{S(T)}. 



Proof. First we simplify 



(19) 

Note that 



E 

k^i 



1 



= lim 



{Xi-XkY 



y 

^ (z- 



T{X)=x 



xy 



(z - Xi 



E 

T(A)=: 



T'{z) 



{z — A) T{z) — X 



□ 



Then 
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That is 



V 1 _ T'\z)-T"{zmz)^x) 

T{X)=x ^ ' \ \ I I 



So passing in a usual way to the hmit in (|19|l we get 
(20) = (A.). 




This means that a diagonal entry of the operator {FD)*{FD) with respect to the 
basis of eigenvectors of J has the form 

□ 

Naturally, in the same way we can find off diagonal entries of the matrix of the 
operator {FD)*{FD). 

Lemma 2.4. For i ^ j 

{{FDYiFD)}, 



T'{K) (A,- A,)2r'(A0 

We would consider A :— {FDF~^)* {FDF^^} as a counterpart of Laplacian due 
to the following proposition. 

Corollary 2.5. A is a positive operator that satisfies 

[J,[J,A]] =2d|0)(0|-2. 
Proof. See Lemma [2.41 □ 

Our plan to estimate [G,J] in ||SJ) was based on the conjecture ||(i^I?)„|| ^ k" 
with K < 1 (typically everything that goes to zero in the subject goes to zero as 
a geometric progression). Since (G„)+ = —{FD)n+ that would give an estimation 
on G: 

||(G„)+||~«"nlogd, 
and we are done. However the following proposition shows that ||(FD)„|| 0. 

Proposition 2.6. There exists the limit 

(21) lim -Ltr{(FD):(FD)J = -^(I - L2)-1L2S(T). 
Proof. Let us use the Chain Rule for the Schwarz derivative 

S{Tn+i) = S{T,,) oTT'^ + S{T). 

Since L^lg o TT'^} = g holds for every function g, we have L^+^{5(r„) o TT'^} = 
L^S{Tn) and therefore 

(22) L^+^S{T^+,) = £^5(r„) + L^+^S{T) = L^SiT) + ... + Vi+^S{T). 

The spectral radius of L2 less than 1 /(P (see Lemma I3.2|l . So (|^ completes the 
proof. □ 
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Remark. We still believe in the limit periodic property of J(/i). Recall that we 
have to estimate not {FD)n itself but the commutator [G„, J„]. Probably it worth 
to mention that the right hand side of the commutant identity for {FD)n, 

{FD)nJn - JniFD)n = F„ - d"|0)(0|F„ 

goes to zero in norm (it's again Lemma [3.2(1 . That is asymptotically {FD)n and 
J„ commute. 

□ 

3. Partial result in the right direction 
3.1. Renormalization equation. Let 

Ty=x 

be a Ruelle operator associated with an expanding polynomial T{z). If J is the 
Jacobi matrix associated with a measure a supported on E, J := J((t), then the 
Renormalization Equation 

(23) V*{z~ jy^V ^{T{z)- J)-^T'{z)/d, V\k) = \kd), 

has a unique solution J := J((t), where a := L*{ct) P], |H|- It follows basically from 
the identity 



1 



z-y 



-{9oT){y)\{x 



T'{z)/d 
T(z) — X 



and the functional representations of both operators in and L^^ respectively. 
Note that becomes JT)) if tr = /z, since for the balanced measure we have 

^i = L*{^l). 

Lemma 3.1. Let J^^'^ be the s-th d x d block of the matrix J , that is 



(24) 



qsd Psd+l 
Psd+1 qsd+1 Psd+2 



Psd+d-2 qsd+d-2 Psd+d-1 
Psd+d-1 qsd+d-1 



Then its resolvent function is of the form 
(25) (o\{z - J^'^)-^\q^ 



T'{z)/d 



T(>>){z) • 

Moreover at the critical points {c : T'(c) — 0} the following decomposition in a 
continued fraction holds true 



(26) 



r(^)(c) = r(c)-g. 



T{c) - qs-i - 

Proof. We write J as a d x d block matrix (each block is of infinite size): 

Qq Vi s+rd^ 



(27) 



J 



Vi Qi 



VdS* 



V2 

'Pd-2 Qd-2 Vd-l 

Vd-i Qd-i 
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Here Vk (respectively Qk) is a diagonal matrix Vk = diaq{pfe+sci}s>o and 5*+ is the 
one-sided shift. In this case V* is the projection on the first block-component. 
Using this representation and being well known identity for block matrices 



we get 



(28) 



T{z)-J 



T'{z)/d 



'A 


B 


-1 r 


C 


D 




z — 


Qo 





{A 



BD- 
* 



'cy 



S+Vd\ {2 

where Ji is the matrix that we obtain from J by deleting the first block-row and 
the first block-column in H27I) . Note that in {z— Ji) each block is a diagonal matrix 
that's why we can easily get an inverse matrix in terms of orthogonal polynomials. 

Let us introduce the following notations: everything related to J*^*^ has super- 
script s. For instance: p^^'' = Psd+k, ^ < k < d, respectively P^'^^ and Q^^^ mean 
orthonormal polynomials of the first and second kind. In this terms equation H28() 
is equivalent to the two series of scalar relations corresponding to the diagonal and 
off diagonal entries 

(29) 



T{z) ~ q^+i 
T'{z)/d 



P, 



\z) 



2 Qd-li^)/Pds 
Pds 



and 
(30) 



Ps+l 



Pi ■ 



■Pd 



T'{z)ld 

We have to remind (see (I25II and (I30|) ) that 



,d-l 



T'{z)/d 



P'd\^) 



(31) 



Now, due to the Wronskian identity, if T'(c) = then 

1 



-PdsQdliic) 



So, combining (j^ . (|Sn|) and we get the recurrence relation 
(32) r(c)-g,+i=T(^+i)(c) + 

with initial data 

r(o)(c) = r(c)-go. 



□ 



3.2. psd„ are exponentially small. 

Lemma 3.2. Let J he the Jacobi matrix associated with iterations {T!„}„>i of an 
expanding polynomial T , that is J — J{fJ.) where L*ii — fi. Then 

(33) C- (pdTps < psd„ < C+ {pdY'ps 

with C± > and < p <l/d. 
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Proof. We recall that Psd = Ps o T and Q^d = {T'/d)Qs o T 0, 0. We use an 
interpolation formula 

(34) ( Rdn= Riy)j^iy)' degR<sd, 
and the Wronskian identity 

(35) Psd{Psd-lQsd - Qsd-lPsd] = 1- 

Substituting (ESJ in ^ we obtain 

Pld^ i {PsdPsd-iY dp= ^ {p,dPsd-i{y)Y'^{y) = X! in \( V 

J . /_ > „ ^sd Wsd^sd)(y) 



Therefore, 



y-Psd{y)=0 y:P,d{y)=0 

E E 
E 



.:P.(.)=0.:T(,)=.((^'Vrf)(QsPi)oT)(y) 



(gsP')(2^) 1 ^ T'2(v) 

£c:P,(a;)=0 " |_ y:T{y)=x ' 

Now we use the Ruelle version of the Perron-Frobenius theorem [7] , |S] with respect 
to L2 (|18|) . According to this theorem 

-^L'^g h{x) j gdv, 

uniformly on x with a certain continuous function ft, > and positive measure i^; 

is the spectral radius of i2- Combining this with the interpolation formula we 
get both-sided estimate (|33|) . 

We only have to show that {pd) < 1. Let b{z) be the complex Green's function 
of the domain C \ E with respect to infinity. Consider the sequence of functions 
{/n}n>i, where fn{z) ■= {b'^"~^Pd„-i){z). It is a multiple-valued function in the 
domain C\E with a single- valued modulus which has a harmonic majorant Un{z): 
|/n(2)P < u„(2;). Moreover, Un{oo) = ||Pd„-i||^2 = 1- We claim that /„ should 

go to zero pointwise. If not then we can find a subsequence {fnt} that converges to 
a non trivial function /. However, in this case, {hf){z) is a non trivial single valued 
in C \ -E function, |(6/)(z)p has a harmonic majorant and (6/) (00) = 0. This 
contradicts to the well-known fact that analytic capacity (that is the Lebesgue 
measure in this case) of E is zero. 

Therefore the sequence converges to zero. In particular 

(6'*"-iPd„_i)(oo) = i -^^0, n^oo. 

Pi---Pd„-i Pi 

But ^ - {pdy\ thus (pd) < 1. □ 
Remark 3.3. Let us mention here that qsd = Qo since 



Qsd 



J yP], dp^ J yP! o TdL*p - J {Ly)P^ dp 
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and Ly = qq. 

3.3. The result. First we prove (undoubtedly well-known and simple) 

Lemma 3.4. Assume that two measures a and a are mutually absolutely continu- 
ous. Moreover, da = f da and 1 — e < f < (1 — e)~^ . Let us associate with these 
measures Jacobi matrices J — J (a), J — J (a). Then for their coefficients we have 

\~ I e . . . . 
\Ps-Ps\ < 



1 - e' 

Proof. Let us use an extreme property of orthogonal polynomials, 
p1..pI = I p1...pIP! da >il-e) I {z^ + ...f da 
>(l-e) inf I P^da={l-e)pl...pl. 



Similarly 



Therefore 



and hence 



2 2 \ /I \ ~2 ~2 
Pl-Ps-1 > (1 - <^)Pl-Ps-l- 

1 



TPi > Pi > (1 - ^ypi 



e 

-eps <Ps-Ps < Ps- 

1 — e 



□ 



Now, we are in position to prove Theorem 1 1.21 
Proof. As it follows from Lemma l3. II 

J T[c) - X 

Here z/*^") is a discrete measure such that suppli/^*-*} C [— ^,^], i^^^'' ?]) — li 
recall that [— ^, ^] is the smallest interval containing the Julia set E. In particular, 

T<->(c)^r(c)-,„-A z^^;^^ 

J T{c) - X 

Now, since 

dist{e:T'(c)=0}{T(c), K,^]} = 5 > 0, 

for every e > there exists n such that 

T[c) - qo 

(here we used Lemma [3 .21) . Recall (|25f) . so. Lemma [3.41 completes the proof. 

□ 
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